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Abstract. Under the assumption that the gas mass frac- 
tion of galaxy clusters estimated out to an outer hydro- 
static radius is constant, it is possible to constrain the 
cosmological parameters by using the angular diameter 
distance relation with redshift. We applied this to a sam- 
ple of galaxy clusters from redshifts of 0.1 to 0.94, for 
which published gas and total masses are available from 
X-ray data. After scaling the gas fraction values to the 
rsoo radius (Evrard, Metzler, & Navarro 1996), we find an 
apparent decrease in gas fractions at high redshifts, which 
can be scaled back to the mean gas fraction value at z ~ 
0.1 to 0.2 of (0.060 ± 0.002) h~ 3 / 2 , when Cl m = 0.55±g;|l 
(Q m + f^A = 1) 1-f statistical error). However, various 
sources of systematic errors can contribute to the change 
in gas mass fraction from one cluster to another, and we 
discuss such potential problems in this method. 
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1. Introduction 

In recent years, measurements of gas mass fraction in 
galaxy clusters together with the universal value for the 
cosmological baryon density as derived from nucleosyn- 
thesis arguments have been used to constrain the cosmo- 
logical mass density of the universe. There are two ba- 
sic assumptions used in such an analysis: (1) that the 
galaxy clusters are the largest virialized systems in the 
universe and based on hierarchical clustering models that 
clusters represent composition of the universe as whole, 
and (2) that the gas mass fraction when measured out 
to a standard (hydrostatic) radius is constant. Evrard 
(1997) applied this argument to a sample of galaxy clus- 
ters studied by David, Jones & Forman (1995) and White 
& Fabian (1995), with redshifts up to 0.2. The gas mass 
fraction values were scaled to the r 50 o radius first defined 
by Evrard, Metzler & Navarro (1996), which has shown 



to be a conservative estimate of the outer hydrostatic ra- 
dius based on numerical simulations. When both gas and 
total (virial) masses are estimated out to this radius, then 
the assumption that all clusters should have the same gas 
mass fraction is well justified. The mean gas mass fraction, 
as derived by Evrard (1997), for the low redshift clusters 
is j^7 ray (r 5 oo) = (0.060 ± 0.003) /i" 3 / 2 . Recently, Myers 
et al. (1997) observed the Sunyaev-Zel'dovich (SZ) effect 
(Sunyaev & Zel'dovich 1980) in a sample of low redshift 
clusters, and derived a mean SZ gas mass fraction, f^ s , 
that ranges from (0.061 ± 0.011) /i" 1 to (0.087 ± 0.030) 
h^ 1 (we refer the reader to Birkinshaw 1998 for a recent 
review on the SZ effect). Using the mean cluster gas frac- 
tion and the nucleosynthesis derived value for the fif,, both 
Evrard (1997) and Myers et al. (1997) put constraints on 
the cosmological mass density of the universe, f2 m , with 
some dependence on the Hubble constant. 

If cluster gas fractions are indeed constant in a sam- 
ple of galaxy clusters when calculated out to the outer 
hydrostatic radius, then it is possible to constrain the cos- 
mological parameters based on their dependence in the 
angular diameter distance relation with redshift. In Sec- 
tion 2 of this paper, we present the potential possibility of 
using cluster gas fraction as a standard candle and apply 
this to a sample of ~ 53 clusters that range in redshifts 
from 0.1 to 0.94, for which published data on masses are 
available from literature. In Section 3, we discuss various 
systematics uncertainties and possible selection effects in 
this method. 

2. Method & Results 

In general, the gas mass of a galaxy cluster is calculated by 
integrating the electron number density profile, n G (r), over 
an assumed shape. The X-ray emission data are usually 
fitted to an apriori assumed distribution, and presently 
the /3-model is used as the preferred option. The gas mass 
within the galaxy cluster is given by: 



M gas (< R)=m 



p I ^c0 
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where n e o is the central number density, r c is the core 
radius of the gas distribution, and (3 is a parameter that 
describes the shape of the distribution (e.g., Cavaliere & 
Fusco-Femiano 1976). For the purpose of this paper, when 
n c0 is determined through the X-ray emission: 



M^(<9)^D A /2 , 



(2) 



while when the SZ decrement (or increment) towards a 
cluster is used to determine n c0 : 



M g s a z s (< 9) cx Dl, 



(3) 



where Da is the angular diameter distance to the cluster. 

The total virial mass of a galaxy cluster can be known 
through the measurements of the intercluster electron 
temperature profile, T c (r), based on the X-ray spectral 
data. Under the assumption of virial equilibrium and as- 
suming an isothermal temperature distribution with elec- 
tron temperature of T e o, the total mass is: 



M tot (<R) 
and 



k B T c0 r rflog n c (r) 
G\xm v d log r 



M tot (< 9)kD a . 



(4) 



(5) 



Apart from X-ray temperature measurements, the total 
mass can be estimated using strong gravitational lens- 
ing, which allows mass measurements near the core re- 
gion of the cluster where background sources are lensed 
to arcs, weak lensing and cluster internal velocity disper- 
sions, which allow mass measurements to outer regions of 
the cluster. However, contrary to X-ray and optical virial 
analysis derived masses, both SZ and lensing masses mea- 
sure the mass along the cylindrical line of sight across the 
cluster. The total masses derived from X-ray temperature 
and strong lensing seem to be off by factors of 2 to 3, while 
weak lensing, velocity dispersion, and X-ray temperature 
masses seem to agree with each other at very large radii. 
The reasons for difference between X-ray temperature and 
strong lensing masses have been studied based on the pres- 
ence of magnetic fields (Loeb & Mao 1994; Miralda-Escude 
& Babul 1995) and cluster cooling flows (Allen 1997). In 
the present paper, we use the X-ray temperature data as 
a measurement of the total (virial) mass, since most of the 
clusters in our sample do not have published lensing and 
optical virial mass measurements yet. 

The two gas fraction measurements based on the X-ray 
total mass is: 



/•X-ray n 3 / 2 
/gas x U K > 

and 

/si « D A , 

for X-ray and SZ gas masses respectively. 



(6) 



(7) 



For Fricdmann-Lemaztre cosmological models, the an- 
gular diameter distance, Da, can be written as: 



D A = 



Hq{1 + z)^/H 



(8) 



where, 



x = ^ [(1 + z'f{\ + fl M z') - z'(2 + z')n A ] 2 dz'.{9) 
Jo 

For Um + Cl\ = 1 (flat universe), \{x) — x an d k = 1, 
for Qm + Cl\ > 1 (closed universe) x( x ) — sin(ar) and 
k = 1 — Qm — Qa, and for Q,m + Q\ < 1 (open universe), 
x(x) = sinh(x) and k = 1 — Om — Oa- When ft a =0, 
the angular diameter distance has a well known analytical 
closed form solution, while for more general case with a 
cosmological constant, the integral over redshift must be 
carried out (Carroll, Press & Turner 1992). 

Using published data, we compiled a sample of galaxy 
clusters and groups with redshifts 0.008 to 0.94. Most of 
the low redshift cluster and group data are from White, 
Jones & Forman (1997) who presented results of an X- 
ray image analysis of 207 clusters based on Einstein data. 
These cluster range in redshift from <~ 0.008 to 0.4, with 
few clusters between z ~ 0.3 and 0.4. Our high redshift 
(z > 0.4) cluster sample is derived from various pub- 
lished data and includes: AXJ2019+1127 (z = 0.94, Hat- 
tori 1997), MS1054-03 (z = 0.829, Donahue et al. 1997), 
MS0451-03 (z = 0.54, Donahue 1996), and C10016+16 
(z = 0.554, Neumann & Bohringer 1996). In order to per- 
form an unbiased test, we calculated both gas and total 
masses in a universe with f2 m = A = and H = 100 
h^ 1 km s _1 Mpc^ 1 , and scaled them to the rsoo radius 
based on relations presented by Evrard (1997): 



rsoo (Too) = (1.24 ±0.09) 



IcO 



10keV 



/i^Mpc. 



At the rsoo radius, the gas fraction is given by: 
/ g a S (r5oo(T o0 )) = /gas(r) ^ r500 | Tc °' ) s 



(10) 



(11) 



where r\ is a scaling parameter determined by numerical 
simulations and r is the radius out to which the gas frac- 
tion is first calculated using mass measurements in lit- 
erature. Evrard (1997) used i] = 0.17, but showed that 
for real data it ranges from 0.13 to 0.17, based on which 
clusters are included. We calculated our gas fraction val- 
ues at r 50 o using r\ values ranging from 0.13 to 0.17, but 
found that the final result presented here is not heavily 
dependent on it. This is primarily due to the fact that for 
massive clusters with cluster gas temperatures close to 10 
keV, the ratio between rsoo and r is within few percent of 
each other. From our total sample of ~ 215 clusters, we re- 
moved all of the clusters with redshifts less than 0.1. This 
is primarily due to the fact that the low redshift sample 
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is mainly made of group size clusters with electron tem- 
peratures less than 3 keV, and the baryonic fraction of 
such clusters are primarily dominated by the baryons in 
galaxies and stars, instead of the intergroup gas compo- 
nent (Fukugita et al. 1997). 

The mean baryonic fraction of galaxy clusters between 
redshifts of 0.1 and 0.15 is (0.060 ± 0.002) h^l 2 , which 
is in agreement with the value derived by Evrard (1997) 
of (0.060 ± 0.003) h~ 3 / 2 . We normalized the cluster gas 
mass fractions to our low redshift mean gas fraction. In 
Fig. 1, we show the observed normalized cluster gas frac- 
tion with redshift. For completeness, we show the whole 
sample, including clusters and groups with redshifts less 
than 0.1. The apparent rise in the gas fraction with red- 
shift at the low redshifts (z <~ 0.008 to 0.1) is primarily 
due to a selection bias: as you go higher in redshift more 
massive clusters are selected in X-ray studies which have 
an increasing gas mass component. 



Fig. 1. The observed gas mass fraction of galaxy clus- 
ters calculated with an angular diameter distance assuming 
Q m — Qa = 0, and normalized the low redshift mean gas 
mass fraction of (0.060 ± 0.002) /*T 3/2 . If the universe is in- 
deed Q m — S!a = then the observed gas fraction throughout 
the redshift plane is expected to be constant. The apparent 
decrease in the gas mass fraction at high redshifts can be ex- 
plained with a cosmological model where O m is low. The drawn 
curve corresponds to Q. m — 0.55 (fl m + Oa = 1)- 



As shown in Fig. 1, the high redshift clusters have gas 
fractions lower than the low redshift sample. The decrease 
in gas fraction at high redshift can be explained in a uni- 
verse with a low Q m . Considering a y 2 minimization to 
the observed data with the angular diameter distance, we 
derived a best fit Q m of 0. 551^23 (l _cr statistical error), 
with Q, m + J1a = 1, at a minimum \ 2 of 62.5 for 53 clus- 
ters. The data are also consistent with Q m < 0.72 (90% 
confidence interval) for an open universe with f^A = 0. 

3. Discussion 

The primary assumption in constraining the cosmological 
parameters using the angular diameter distance is that 



the gas fraction at high redshift is same as the fraction 
at low redshift clusters. It is likely that gas fractions may 
vary from one cluster to another, but numerical simula- 
tions seem to strongly suggest that gas fractions measured 
to the outer hydrostatic radius is constant and is not 
significantly influenced by hydrodynamical effects (e.g., 
Evrard, Metzler & Navarro 1996; Evrard 1997; Lubin et 
al. 1996). In Shimasaku (1997), the assumption of con- 
stant gas (baryonic) fraction was used to put constrains 
on cr 8 , the rms linear fluctuations on scales of 8 h _1 Mpc, 
and on n, the slope of the fluctuation spectrum. However, 
possible reasons for fluctuations in the gas fraction needs 
to be considered, and we discuss briefly various systematic 
uncertainties and biases in this method. These include gas 
clumping, presence of magnetic fields and magnetic pres- 
sure support, various non-gravitational processes such as 
heating and cooling in galaxy clusters, additional baryonic 
contributions, and a possible phase of gas injection from 
galaxies in clusters to intercluster medium. 

3.1. Gas Clumping & Magnetic Fields 

Given that the gas masses are estimated from X-ray emis- 
sion, which measures the < n 2 >, one needs to assume 
that the gas is smoothly distributed within the interclus- 
ter space. However, if gas is clumped with a clumping 
factor C (=< n 2 > x l 2 / < n e >), then the amount of 
gas required to produce a fixed X-ray emission will re- 
duce by C. The observational evidence against clumping 
is presented in Evrard (1997). Unlike the X-ray gas mass, 
the gas mass derived using the SZ effect is not subject to 
gas clumping as it measures < n e > directly (although 
the Hubble constant that could be derived by combining 
SZ and X-ray data is subject to gas clumping effects). 
Therefore, one should be able to constrain clumping in 
galaxy clusters based on a set of SZ and X-ray data for 
a cluster sample. Currently, there are 4 mass estimates 
from SZ effect (Myers et al. 1997), with a mean gas frac- 
tion of (0.087 ± 0.030) hr 1 to (0.061 ± 0.011) hr 1 . If 
we used the mean X-ray derived gas fraction, and un- 
der the assumption that the difference between X-ray and 
SZ gas fractions is due to a combination of gas clumping 
and the Hubble constant, an upper limit to the clump- 
ing factor C is <~ 1.24 when h — 0.65. There is also the 
possibility that gas clumping is nonisothermal such that 
C =< n\T e > 1 / 2 / < n e >< T e >, but, constraints on 
such a gas clumping model will only be possible after re- 
liable cluster gas temperature profiles are measured. In 
future, this should be possible with planned X-ray mis- 
sions such as AXAF, XMM and Astro-E. 

The presence of magnetic pressure support has been 
discussed as a solution to explain the difference between 
X-ray virial and strong gravitational lensing mass in the 
core regions of the cluster. Given that our method uses the 
fact that gas fraction out to an outer hydrostatic region is 
constant, it is unlikely that magnetic fields play a major 
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role in varying the cluster gas fraction from one cluster to 
another. However, the difference in gas fraction near the 
core regions of the cluster may be explained using mag- 
netic fields. Here again, present data do not allow reliable 
estimates on the amount of magnetic pressure support in 
galaxy clusters. 

3.2. Heating & cooling in galaxy clusters 

Cluster cooling flows have been observed up to ~ 40% of 
the X-ray clusters, and it is expected that ~ 10% of the 
observed X-ray emission in such clusters are due to cool- 
ing flows (Fabian 1994). Such cooling flows are expected 
to increase the gas fraction near the core region of the 
cluster where such effects exist. However, the apparent in- 
crease in the cluster gas fraction from high to low redshift 
cannot be fully explained as due to cooling flows, since the 
low redshift cluster sample both contains high cooling flow 
and low cooling flow clusters, with gas fractions that agree 
with each other, and that there is little evidence for cool- 
ing flows in the high redshift sample. If a heating process 
exists within galaxy clusters then the cluster gas fraction 
may be lowered. However, no physical mechanism has yet 
been suggested which can lower the cluster gas fraction. 

3.3. Additional baryonic contributions 

The arguments used here completely ignored the presence 
of baryons other then gas in clusters. It is well known that 
the baryonic component of the massive clusters are dom- 
inated by the X-ray emitting cluster gas rather than the 
baryons associated with stars and galaxies in clusters. Ac- 
cording to White et al. (1993), the mass ratio of stellar 
component to the X-ray emitting gas is 0.2 hr 1 for the 
Coma cluster. Assuming Coma is a typical rich cluster, 
we can assume that all other clusters have the same ratio. 
However, Hattori (1997) suggested that AXJ2019+1127 
at a redshift of 0.94 is a dark cluster with mass to light 
ratio M/L ~ 1500 h Mq/Lq. Therefore, the stellar com- 
ponent may not be constant from one cluster to another, 
but, the effect would be to decrease the baryonic fraction 
at low redshifts from the high redshift values, contrary to 
what is observed in the present data. Also, various forms 
of baryonic components are expected to exist in the inter- 
cluster medium (e.g., intercluster stars, Uson et al. 1991), 
but, such contribution is expected to be lower than the 
observed variation in the gas fraction. 

4. Conclusions 

Under the assumption that the baryonic (gas) fraction in 
galaxy clusters are constant, we have constrained the cos- 
mological parameters using the angular diameter distance 
relation with redshift. However, the present X-ray based 
data from various studies may be providing a biased re- 
sult due to unknown selection effects. When the angular 
diameter distance dependence on the traditional Hubble 



constant measurements using SZ and X-ray data are in- 
cluded with the present method involving cluster gas mass 
fraction, it is likely that tighter constraints on the cosmo- 
logical parameters are possible. We hope to explore this 
possibility in a future paper. Given the biases that may 
go in to the presented gas mass fraction vs. redshift di- 
agram, we suggest that results from a well defined and 
random sample of galaxy clusters, such as that would be 
available from AXAF in X-ray and from intcrferometric 
observations in SZ, be used to constrain the cosmological 
parameters using the angular diameter distance relation 
with redshift. 

5. Note 

Since this work was first submitted, a paper by Danos 
& Pen (1998) appeared, arriving at qualitatively similar 
conclusions - though using cluster electron temperature, 
luminosity and angular size to measure the gas mass frac- 
tion of three high redshift clusters, which were compared 
to a global gas mass fraction value determined from nu- 
merical simulations. 
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